Transport Through Nanostructures with Asymmetric Couphng to the Leads 
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Using an approach to open quantum systems based on the effective non-Hermitian Hamiltonian, 
we fully describe transport properties for a paradigmatic model of a coherent quantum transmitter: 
a finite sequence of square potential barriers. We consider the general case of asymmetric exter- 
nal barriers and variable coupling strength to the environment. We demonstrate that transport 
properties are very sensitive to the degree of opening of the system and determine the parameters 
for maximum transmission at any given degree of asymmetry. Analyzing the complex eigenvalues 
of the non-Hermitian Hamiltonian, we show a double transition to a super-radiant regime where 
the transport properties and the structure of resonances undergo a strong change. We extend our 
analysis to the presence of disorder and to higher dimensions. 

PACS numbers: 05.50.-|-q, 75.10.Hk, 75.10.Pq 



I. INTRODUCTION 

Open quantum systems, which exchange matter and 
energy with an environment, are at the center of many 
research areas in condensed matter, atomic, molecular, 
and nuclear physics. Major problems of current interest 
range from quantum computing to transport in meso- 
scopic systems to basic theoretical issues, including the 
measurement problem in quantum mechanics. 

The nature and degree of opening affects the properties 
of a system in a highly nontrivial manner. An example of 
this is the super-radiance phenomenon in a finite quan- 
tum system coupled to an environment characterized by 
a continuum of states. Generically, at weak coupling, all 
internal states are similarly affected by the opening and 
acquire small decay widths, resulting in narrow transmis- 
sion resonances. As the coupling increases and reaches 
a critical value, the resonances overlap, and a sharp re- 
structuring of the system occurs. Beyond this critical 
value, a few resonances become short-lived states, leav- 
ing all other (long-lived) states effectively decoupled from 
the environment. This general phenomenon is referred to 
as the super-radiance transition [l|, Q , due to its analogy 
with Dicke super-radiance [3] in quantum optics. 

In a recent work [ij generalizing the schematic tight- 
binding model discussed in 0, Q, it was shown that 
the phenomenon of super-radiance actually occurs in the 
problem of transport through realistic nanosystems. The 
specific situation analyzed in Ref. jj] is transport through 
a one-dimensional (Id) sequence of potential barriers, see 
Fig.[TJ This paradigmatic model of solid state physics ap- 
pears in many important applications, including semicon- 
ductor superlattices and one-dimensional quantum dot 
arrays. It has been widely discussed in the literature [Q- 



|8|; the transport properties have been analyzed as the 
system-environment coupling was varied by adjusting 
the widths of the external barriers. In Ref. j4i| only sym- 
metric coupling was considered, i.e. equal left and right 
external barriers. It was shown that maximum transmis- 
sion through the array is reached precisely at the super- 
radiant transition. The transport properties of a Id se- 
quence of potential barriers were analyzed with the aid of 
the energy-independent effective non-Hermitian Hamilto- 
nian. This approach produces excellent agreement with 
the exact (numerical) treatment of the problem for weak 
tunneling between the wells. 

In this paper we use the same framework to analyze 
the transport properties of a Id sequence of potential 
barriers with asymmetric coupling to the environment. 
We vary the external coupling strength, keeping the ratio 
between left and right coupling constant. This allows us 
to determine the maximum transmission as a function 
of the asymmetry and to extend the analysis to higher- 
dimensional systems: quasi-ld, 2d, and Sd. 

The paper is organized as follows. In Sec. [Ill we intro- 
duce the energy-independent effective Hamiltonian that 
relates the sequence of barriers to the open Anderson 
model, and show how the transmission properties can be 
determined in this formulation. In Sec. IIIIl we find the 
strength of coupling to the environment at which both 
integrated transmission and average transmission at the 
center of the energy band are maximized, and derive the 
scaling of this critical coupling strength as a function 
of asymmetry and degree of disorder. The structure of 
resonances is analytically computed as a function of en- 
ergy and of asymmetry of the coupling for the case of no 
disorder. Analyzing the complex eigenvalues of the non- 
Hermitian Hamiltonian in Sec. IIVI a double transition to 



a super-radiant regime is shown to occur. Contrary to 
the symmetric coupling case, the maximum transmission 
in the asymmetric case is not reached at either transition, 
but occurs instead at a critical value located between the 
two. The super-radiant transitions have significant con- 
sequences for the transport properties of the Id Ander- 
son model. The number of resonances decreases by one 
every time a super-radiant transition is crossed. For a 
large number of sites in the Id chain, a clear signature 
of the two super-radiant transitions is observed when we 
analyze the resonance structure near the center of the 
energy band as a function of the coupling strength to the 
leads. In Sec. |Vl we compare our results with the ran- 
dom matrix theory of transport, showing that random 
matrix theory is only partially applicable to the Id An- 
derson model. Finally, in Sec. IVIl we extend our results 
to multi-dimensional scenarios, showing that the validity 
of our findings is not limited to the Id case. 



II. EFFECTIVE HAMILTONIAN 

The effective Hamiltonian approach to open quantum 
systems was formulated in the book [9] for nuclear re- 
actions, and later generalized and studied in detail, see 
for example [l!, St i, i, [13, [HI. In Ref. j^] two of 
the present authors demonstrated how to build an effec- 
tive non-Hermitian Hamiltonian that correctly describes 
transport through a sequence of potential barriers. Here 
we will only review the main points of the approach and 
establish definitions and notations. 

In the absence of disorder, we consider quantum trans- 
port through a sequence of A^ -I- 1 potential barriers of 
height Vq and inter-barrier separation L, as illustrated in 
Fig. [H All A^ — 1 internal barriers have width A, while 
the two external barriers have widths Ai.2. 
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Sequence of potential barriers of finite height and 



We may compute transmission through the system in 
a standard way, by matching the wave function and 
its derivative on either side of each barrier in Fig. [T] 
Throughout the paper, we will use units with h jlrrif, — 
1. Thus, when distances A, Ai_2, and L are measured in 
nm (the typical scale for semiconductor superlattices), all 
energies are calculated in units of 38 meV. In what follows 
we set L = 2, A = 0.2, Vb = 1000, and E^ = VqI2 = 500 
(so that the energy shift vanishes, see below). 

In the limit of weak tunneling between the sites, a se- 
quence of A^ potential wells behaves as an open Id Ander- 
son model. As shown in Ref. [4J, the Id effective Hamilto- 
nian for an energy band centered at Ei^ can be written in 
the site basis in a way similar to that used in Refs. 2, 5], 
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The edge states, |1), localized in the first well on the left, 
and |A^), localized in the last well on the right, acquire 
finite widths, 712, and energy shifts, 61^2, due to the 
coupling to the environment. By inter-site tunneling, this 
coupling propagates through the chain. 

The effective Hamiltonian Ji correctly reproduces 
transmission through a sequence of potential barriers if 
we define the tunneling coupling n as 



n = 



2a^Eo 



K)(l + aL/2) 



exp (— aA) , 



(2) 



where a = ^Vq — Eq. Similarly, the widths and energy 
shifts can be written as 



71.2 = T#Tf^ exp (-2aAi,2), 



(3) 



<5l,2 = 



iak 



■71,2: 



where k = ^/Eq. The shifts ^1^2 vanish for Eq — Vb/2; 
otherwise the sign of S is given by the sign oi Eo — Vo/2. 
We will study how the transport properties depend on the 
system-environment couplings 71,2, which are varied by 
adjusting the external barrier widths A1.2 while keeping 
all other parameters fixed. 

With the aid of the effective Hamiltonian, the trans- 
mission coefficient T°'^{E) from channel b to channel a 
can be determined. 



T''\E) = \Z''\E)[' 



where 
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is the transmission amplitude. The channels are labeled 
by the quantum numbers that characterize the contin- 
uum states in the environment, not including the energy. 



In the Id case we have two channels, a = 1,2, correspond- 
ing to the left and right scattering states, respectively. 
The factors A'^ represent the transition amplitudes from 
state \i) to channel a. In our arrangement, the only 
non- vanishing transition amplitudes are ^} = ^/ji and 
Aj^ = y/j2- The complex eigenvalues £k of H coincide 
with the poles of Z{E). The spectrum of the complex 
eigenvalues of the effective Hamiltonian is of great im- 
portance for understanding the transport properties of 
the system. 

Using the effective Hamiltonian of Eq. ([T|), the trans- 
mission between left and right leads can be computed 
as 



by 



T{E) 



(V7i72/f^) 
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m=iiE-£k)/^ 



(6) 



where here and in the following T{E) = T^'^{E). It was 
shown in Ref. [4] that the exact transmission obtained 
by matching of the wave functions is in excellent agree- 
ment with the effective Hamiltonian approach, Eq. ([6]), 
for aA > 1. 



III. TRANSMISSION THROUGH 

NANOSTRUCTURES WITH ASYMMETRIC 

COUPLING 

In this Section we analyze the behavior of the trans- 
mission as we increase the coupling to the external envi- 
ronment, keeping the ratio of the two external couplings, 
7i = 7 and 72 = 7/9, fixed and equal to q, the asymmetry 
parameter. We treat first the case of no disorder, where 
the unperturbed states in all sites have the same energy 
Eq, and later extend to the disordered case. We will also 
drop the energy offset Eq from the effective Hamiltonian, 
so that the center of the energy band is always at zero 
energy, and for simplicity we neglect the energy shifts ^1^2 
in the following considerations (which in the potential of 
Fig. [1] corresponds to Eq = V/2). 



A. Ordered Case 

Before considering the problem of a general A^-level 
system, we first treat the one-well and two-well cases, 
which give us valuable insight into how transport prop- 
erties change as we increase the coupling strength to the 
environment. 

We start with the transmission through a single quan- 
tum level, namely through a quasistationary state cre- 
ated by two potential barriers. In this case we have only 
one resonance and the effective Hamiltonian reduces to 
"H — —i'y/2 — i'y/2q. The resonance height is independent 
of the coupling 7, and transmission is never perfect for 
asymmetric barriers. Indeed, from Eq. (jS]) we see that 
maximum transmission is attained at E = and is given 



T{E = 0) 



4q 



{q+iy 



(7) 



so that T{E = 0) = 1 only for the equal-coupling case of 
q=l. 

The situation is different when we consider transmis- 
sion through two quantum states. This problem has been 
studied previously, see for instance Refs. [12], [ij] and ref- 
erences therein. In fact, many of the results obtained for 
the two-level case are of more general validity. The effec- 
tive Hamiltonian for two originally degenerate levels can 
be written as 



n = 



-11 



-f7/g 



(8) 



From Eq. ^ we can then compute the transmission. In 
particular, the transmission at _B = is 



T{E ^ 0) 
At the critical value. 



h/^)/V^ 



1 + {i/n)y4q 



O -^^' 



(9) 



(10) 



transmission is perfect. This is in contrast with the one- 
level situation, where transmission is never perfect for 
q^l. 

By applying the residue method to Eq. ^ , we can also 
compute the normalized integrated transmission; 



S 



1 

in 



T{E) dE 



TT"f/Q 



2(g+l)[l + (7/r!)V4<z] 



(11) 



where 4f2 is the width of the energy band. The maxi- 
mum integrated transmission occurs at the same critical 
7 value given by Eq. (|10p . The quantitative behavior of 
S is important in applications, for instance in the design 
of electron band-pass filters for semiconductor superlat- 
tices [6|. 

We now analyze transport properties for the general 
N -level system, many of which parallel those of the spe- 
cial one- and two-level cases. Numerically calculating 
transmission as a function of energy, we find that perfect 
transmission is attained at all resonance peaks precisely 
at the critical coupling given by Eq. (|10p independently 
of the value of N, see for example the upper right panel of 
Fig. [6] Moreover, for all A^ > 2, the integrated transmis- 
sion has a maximum at the same value of the coupling. 
The simple theoretical expression pTj) . obtained for the 
integrated transmission in the two-level case, closely re- 
produces the integrated transmission for any N > 2, see 
Fig. [2] The origin of this iV-independent behavior in the 
coherent transport regime has been discussed previously 
in the context of symmetric coupling [J| . 
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FIG. 2: (Color online) Integrated transmission S obtained 
in the effective Hamiltonian approach for A'' = 100 sites and 
different asymmetry parameters q: squares refer to g = 1, 
circles to g = 10, and diamonds to g = 25. The solid curves 
represent the result (|lip obtained for N = 2. The vertical 
dashed lines indicate the critical value of the coupling strength 
in each case, "f/^l = 2y^, where 5* is predicted to have its 
maximum, see Eq. (|10p . 



Our results agree with the full analytical expression 
for the transmission amplitude ([5]) in the general A^-level 
case [2|, 



Z{E) = - 



where 



2^jT^P-{E) 



1 + ^(7l + 72)P+{E) + 7172[P2 (S) - Pl(E)] ' 



(12) 



1 ^ 1 

p^(E) = — y (±1)" 

^ ' A^ + iZ^^ ' E-En 



sni 



iV+1 



(13) 
and the sums are over the unperturbed Bloch wave ener- 
gies of the closed system. 



En = 2n cos 



N + 1 
This determines the exact transmission probability, 

471 72 -P^ 



T{E) 



l+lll2[Pl-Pl)\ +{11+12? Pi 



(14) 



(15) 



Inside the energy band, at any pole E — En corre- 
sponding to a Bloch eigenstate, both sums P± diverge 
with |P+/P_| -^ 1, and the transmission takes the A^- 
independent and 7-independent value given by Eq. ([7]), 



T{E - En) - Ti 



4g 



(9 + 1)= 



(16) 



obtained above for £■ = in the special case N = 1. 
Outside the band, for \E\ > 217, the sum P+ converges 
at large N to the ^''-independent value 



P+{^) 



2n 



(17) 



(see Appendix), whereas the sum P 
small at large N, 

P-(e)«T 

where e — E/2VI and z± = 
cay of P- results in exponentially weak transmission. In 
Eqs. (fT7|) and ([TS]). the upper and lower signs should be 
chosen for e > 1 and e < — 1, respectively. For energies 
E inside the band, the sums in Eq. (IT51) take a simple 
form, see Appendix: 
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f sin(JV/3) 
p,IA ~ _ J sin[(W+l),3] 

I sin[(Ar+l)/3] 



(19) 



where /3 = cos~^(e). 

A convenient simplification does occur for energies near 
the middle of the band, \E\ <^ il, where for large N 
the sums P± are dominated by terms associated with 
n w A'^/2, while distant contributions from the left and 
right sides of the Bloch spectrum cancel. Near the center 
of the band. En may be replaced by a picket fence spec- 
trum with spacing D — 2tiQ,/N , while the last factor in 
Eq. ([T3| reduces to unity. Defining r by £■ = £'„,, + rD, 
we then have the A^— independent result 



P±{E) 




and therefore 

T{E) 



cot(7rr) 
(±l)"''csc(7rr) 

1 671 72 



(20) 



[AVL + 1112/^Y sin'^(7rr) + 4(71 -|- 72)^ cos2(7rr) 

(21) 
At energies in the Bloch spectrum (r = 0), Eq. ((2T|) re- 
duces, as it must, to the exact expression ([T6|) . while 
midway between any two neighboring poles [r = 1/2) 
the transmission becomes 



T{E = En + D/2) « r2 = 



iim/^q 



1 + (7/r!)V4<Z 



(22) 



which agrees with the transmission at E' = obtained 
above in the special case of A^ = 2 wells, Eq. ([9|). At 
these energy midpoints, perfect transmission occurs at 
the critical value of the coupling given by Eq. (fTO)) . just 
as it does in the special case N = 2. Averaging Eq. (PT|) 
over an energy window containing multiple resonances, 
il/N ^ AE <C ri, we obtain the energy-averaged trans- 
mission near the middle of the band, 



T 



T[E) dr 



81/n 



(<Z + l)(4+(7/f7)2/q) 



(23) 



that reaches its maximum, T™^^ — 2y/q/{l + q), again at 
the critical value of the coupling j/il given by Eq. ([TOl) . 
The critical value associated with the maximum trans- 
mission in Fig. [2] is also consistent with the results of 
Ref. JL4i]. There, the authors found that given a se- 
quence of potential barriers of width A, adding two ex- 
ternal barriers satisfying A^ -I- Ajj — A produces an in- 
crease in transmission while leaving the resonance en- 
ergies unchanged. From Eqs. ([2]) and ([3]), we can see 
that the maximum transmission condition (|10p obtained 
from the tight-binding model coincides with the condi- 
tion A^ -|- Afi — A for the case of E^ = Vo/2, in which 
case the energy shift is zero. 
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B. Disordered Case 



In this subsection we analyze the effect of disordered 
on-site energies. The survival of the super-radiant re- 
structuring in the disordered chain was established in 
Ref. |5]. In Ref. [J] we showed that the effective Hamil- 
tonian, Eq. ([J), correctly describes the sequence of po- 
tential barriers in the presence of disorder (e.g., vari- 
able well width), when the disorder is sufficiently weak. 
We consider random variations of the diagonal energies, 
Eq + SEq, where SEq is uniformly distributed in the in- 
terval [—W/2,+W/2], and W is the disorder parameter. 
In Ref. [3l it was shown that for a A ^ 1, a random 
variation of SEq in the interval [—W/2, +W/2] corre- 
sponds to a random variation of the well width SL in 
[-WL/4:Eo,+WL/4:Eo]. 

The effective non-Hermitian Hamiltonian with diago- 
nal disorder is equivalent to a Id open Anderson tight- 
binding model [ISl . [l6| . The eigenstates of the Anderson 
model are exponentially localized on the system sites, 
with the tails given by exp(— x/Lioc)- Here the localiza- 
tion length Lioc depends on the disorder strength p7l |. 
and X is distance in the direction of transmission, mea- 
sured in units of the well size. For L\oc ^ N, the trans- 
mission decays exponentially with N; this is the localized 
regime. 

Let us first consider the case when the mean level 
spacing D (at the center of the energy band) is not 
strongly modified by the disorder. In the ordered case, 
D « 2tiQ,/N for _E = 0, while for strong disorder the band 
width is W and we have D k, W/N . Thus, we expect that 
for W < 27rJ7 the mean level spacing is not strongly in- 
fluenced by the disorder. This regime is shown in Fig. |31 
where we plot the average transmission as a function of 
the coupling strength 7/fi for N = 100 and W/fl = 0.5. 
As indicated by the vertical dashed lines, the maximum 
transmission is reached at the same critical values of 7/J7 
obtained previously in the absence of disorder, Eq. ([TU]). 
At this critical value, each transmission curve intersects 
the transmission curve for the symmetric case q — 1. In- 
deed, at the critical coupling, the tunneling probabilities 
from the left and right are equal, implying that for this 
value the asymmetric system behaves as a system with 



FIG. 3: (Color online) Average transmission near E — 
as a function of coupling strength 7/n for a disordered Id 
chain of A'^ = 100 wells, for different values of the asymme- 
try parameter q. The disorder strength is W/^ — 0.5. Each 
value of transmission corresponds to the average over an en- 
ergy window —0.1 < E < 0.1, and ensemble average over 100 
realizations of the disorder. The curves are obtained from 
the effective Hamiltonian H (open Anderson model), whereas 
symbols are computed for the sequence of potential barriers. 
Vertical dashed lines indicate the critical values of the cou- 
pling flO)) for the ordered case. 



symmetric coupling, see Sec. |Vl 



The situation is different for strong disorder, W > 
27rf2. The critical value for maximized transmission now 
depends on the disorder strength. For strong disorder we 
enter the localized transport regime studied in Ref. jl8| . 
where transmission is log-normally distributed, so that it 
is more convenient to consider the average of In T rather 
than the average of T. Numerically we find that the 
value at which transmission is maximized is proportional 
to the mean level spacing D at the center of the energy 
band: 7 oc ND, see the arrows that indicate 7 = 1.2ND 
for each curve in Fig. 21 where q = 10. Since D « W/N 
when disorder is strong, the critical value of the coupling 
is in this case proportional to the disorder strength. 

Additionally, we note the close agreement between 
symbols and solid curves in Figs. [3] and 01 demonstrating 
good correspondence between the exact barrier problem. 
Fig. [I] and the energy-independent effective Hamiltonian 
of Eq. ([T]). This correspondence persists for weak and 
strong coupling to the environment, in highly symmetric 
and highly asymmetric situations, and also for both weak 
and strong disorder. 




FIG. 4: (Color online) Average of the logarithm of transmis- 
sion at _E = as a function of the coupling 7/f2 for the Id case, 
with A'' — 100 sites, and asymmetry parameter q = 10. In our 
units, the tunneling coupling is f2 = 1, and the three data sets 
correspond to different disorder strengths W. {lnT)o repre- 
sents the value of the transmission at 7/SI = 0.01, which is 
the lowest value of 7/^ computed. For the sake of clarity, 
we plot (InT) — (InT)o so that all data start from the same 
point. In each case, symbols refer to the results obtained from 
the effective Hamiltonian H: circles for W — i, squares for 
W — 10, and crosses for W = 15, while the curves refer to a 
sequence of potential barriers: dot-dashed for W — 4, dashed 
for W — 10, and solid for W — 15. Averaging is performed 
over 1000 realizations of the disorder. The dashed vertical 
line indicates the critical value of the coupling, Eq. (|10|l . at 
which transmission is maximized for low disorder. The arrows 
indicate 7/fi = 1.2ND/il, where D is the mean level spacing 
at the center of the energy band and the factor 1.2 is obtained 
from fitting. 



IV. DOUBLE SUPER-RADIANT TRANSITION 
AND STRUCTURE OF RESONANCES 

In Ref. [4| we showed that maximuni transmission is 
reached at the super-radiant transition in the symmet- 
ric coupling case. The transition is signaled by a seg- 
regation of resonance decay widths, i.e., of the imagi- 
nary parts of the eigenvalues of the effective Hamiltonian, 
into two groups, super-radiant (short-lived) and trapped 
(long-lived) |i|, 0, Ij, \^- The number of super-radiant 
states is equal to the number of channels coupling the 
system to the environment [two for Fig. [T] and for the 
effective Hamiltonian of Eq. ((l}]. In order to identify the 
super-radiance transition we compute the average of the 
non-super-radiant widths, i.e. of the smallest A^— 2 decay 
widths of the effective Hamiltonian. In Fig. [5l we plot 
this average value as a function of the coupling 'y/il for 
two asymmetry values, q = 4 and q = 10. In each case 
the average width shows two maxima as 7/^2 is varied. 
At each maximum, one of the widths segregates from the 
others. Thus, in the presence of asymmetry, we have 



two super-radiant transitions associated with two criti- 
cal values of the coupling. According to the resonance 
overlap condition, 71,2/^ = 2 ^], the two transitions are 
predicted to occur at the values 

(7/^)sRi = 2, (7/^)sR2 = 2q, (24) 

indicated by the vertical dashed lines in Fig. [SI in very 
good agreement with the numerical results for different 
values of the asymmetry. 




FIG. 5: (Color online) Average width (computed over the 
smallest N — 2 widths) in units of the mean level spacing at 
the center of the energy band. The average width is calcu- 
lated as a function of j/fl for TV — 100 and two values of 
the asymmetry parameter q. Dashed vertical lines indicate 
the locations of the two super-radiant transitions, as given by 
Eq. (dH). 

Comparing Eqs. pO)) and ((24)) . we observe that, in the 
presence of asymmetry, the point of maximum transmis- 
sion is always in between the two super-radiant transi- 
tions. The transmission maximum occurs approximately 
at the minimum of the average decay width, while the 
super-radiant transitions occur at the two maxima of the 
average decay width. This can be compared to the case 
of symmetric barriers, where the maximal transmission 
precisely coincides with the (single) super-radiant tran- 
sition [J], in agreement with Eq. ([M|) for g = 1. 

Since the super-radiant states are very broad, the num- 
ber of observed resonances changes after each transition. 
As we increase the coupling, the number of resonances 
changes from N at weak coupling to A^ — 1 after the first 
transition, and finally to A^ — 2 for strong coupling. This 
change in the resonance structure can have important 
consequences for experimental current- voltage curves, for 
example in semiconductor super lattices [3] . For symmet- 
ric coupling, the number of resonances changes directly 
from A^ to A^ — 2 at the (single) super-radiant transi- 
tion [i|. 

An important consequence of the two super-radiant 
transitions can be seen in the structure of resonances near 
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FIG. 6: (Color online) Transmission for A'^ = 100 and asym- 
metry parameter q = 10 calculated as a function of energy in 
a small interval around the center of the energy band. Ten 
different values of the coupling y/Q, are shown for the case of 
no disorder, W = 0. The horizontal dashed line, Eq. (|16p . is 
the analytically obtained transmission at the location of the 
poles, which is also the maximum transmission for A'' = 1, 
Eq. 0. 



the center of the energy band. In Fig. [5] we consider aid 
chain with N — 100, q — 10, and no disorder {W — 0), 
and calculate transmission as a function of energy for 
different values of the coupling "//^. In each panel in 
Fig. ini the dashed horizontal line indicates the transmis- 
sion at energy values belonging to the Bloch spectrum, 
Eq. (J16p . which is also the maximal transmission attain- 
able for transport through one level only, as given by 
Eq. 0. 

In the upper left panel of Fig. [SJ the resonance struc- 
ture is shown for several cases of weak coupling, 'y/fl < 2, 
i.e. below the first super-radiant transition. The maxi- 
mum resonance height is here independent of 7/^2, and 
equal to the single-level resonance height. Only the 
widths grow as 7/51 is increased in this regime. In the 
upper right panel, we consider 2 < j/^ < 2q, the in- 
termediate regime between the two super-radiant transi- 
tions. Perfect transmission is reached for all resonances 
precisely at the geometric mean of the two transitions, at 
7/r2 = 2y/q. In this regime, it is the minimum transmis- 
sion that is independent of j/ft and equal to the single- 
resonance transmission height, Eq. ^. In the lower left 
panel the region around the second super-radiant transi- 
tion is shown. As •y/il crosses the value j/il = 2q, the 
minimum transmission drops below the one-level trans- 
mission value. Finally, in the lower right panel, we show 
the regime of large coupling, 7/i7 > 2q. Here it is again 
the m,axim,um, transmission that is independent of 7/ri 
and given by Eq. ([71), just as in the weak coupling regime, 
but the resonance widths now shrink with increasing 7/fi. 
The behavior shown in Fig. [HI indicates that, at least for 



large N, a qualitative change in the resonance structure 
near the center of the energy band occurs at each of the 
two super-radiant transitions. 
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FIG. 7: (Color online) The case of iV = 100 sites with 
g = 10 is considered. Upper panel: the resonance height 
(transmission maximum) averaged over the center of the band 
(—0.1 < £ < 0.1), as a function of 7/O. Lower panel: the 
transmission minimum averaged over the same energy win- 
dow. The two vertical lines indicate the positions of the two 
super-radiant transitions. The solid curve shows the predic- 
tions of Eqs. (|25p and (|26|) . A clear change in the transport 
properties is seen at each of the two transitions. The effect 
persists in the presence of disorder {W/Q — 0.2). In each 
panel, the results obtained from the effective Hamiltonian {'H) 
are compared with those obtained from the sequence of po- 
tential barriers (PB). 



To better study this change in the resonance structure, 
we now compute (Tmax) , the average of the transmission 
maxima, and (Tmin) , the average of the transmission min- 
ima, as functions of the coupling strength. In each case, 
the averaging is performed over a window near the cen- 
ter of the energy band, —0.1 < E < 0.1. The upper 
panel of Fig. [7] shows that (Tmax) is constant to the left 
of the first super-radiant transition and to the right of 
the second transition (shown as vertical lines). The con- 
stant value coincides with that for the one-level maximal 
transmission, Eq. ([7]), as shown by the horizontal line. 
Between the two transitions, (Tmax) reaches a maximum 
value at 7/O = 2y/q. In the lower panel of Fig. [71 we 
show the behavior of (T,„i,i) as a function of 7/fi. This 
average rises to reach the one-level transmission value 
at the first transition, stays constant up to the second 
transition, and then decreases. 

For the ordered case we then find numerically the sim- 
ple results: 



\-^ max/ — 



T2 for 2 < -f/n < 2q 
Ti otherwise 



(25) 
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and 



\-^niin/ 



Ti for 2 < 7/f] < 2<7 
T2 otherwise 



(26) 



which agree well with the analysis in Sec. IIII Al Here Ti 
is the transmission at resonant energies in the Bloch spec- 
trum, Eq. (J16p . that coincides with the one- level maximal 
transmission value, Eq. ([7]). Similarly, T2 is the trans- 
mission value at midway points between these energies, 
Eq. dm, that coincides with T{E = 0) for the two-level 
case, Eq. ^. These formulas, represented by the solid 
lines in Fig. [71 are in good agreement with the numerical 
calculation. 

It is also interesting to point out that from the analyti- 
cal expressions, Eqs. (|22)) and (fT6)) . we have that T2 > Ti 
for 2 < 7/n < 2g, so that we regain the thresholds for 
the two super-radiant transitions. 

In the same Figure, we also illustrate the effect of 
adding disorder. For weak disorder, the behavior is simi- 
lar to the ordered case, while for strong disorder the effect 
of the two super-radiant transitions is smoothed out. In 
Fig.[7]the results obtained from the effective Hamiltonian 
(circles and squares) are compared with the results from 
the sequence of potential barriers (crosses and pluses). 
Again, the agreement is excellent, both with and with- 
out disorder. 



V. COMPARISON WITH RANDOM MATRIX 
THEORY 

It is interesting to compare our results on the Id Ander- 
son model with standard results obtained in the frame- 
work of the random matrix theory of transport. In 
Ref. [20| , the dependence of the conductance on the de- 
gree of opening was analyzed for a quasi- Id system with 
symmetric coupling to the leads, while the case of asym- 
metric coupling was later addressed in Ref. [T^ , with dif- 
ferent tunneling probabilities for the right and left ends. 
The tunneling probability used in Ref. [18] corresponds to 
the "elastic scattering" probability r defined in Ref. [21| , 
where the case of symmetric coupling and varying degree 
of internal disorder was considered. This probability for 
channel a is defined as t° = 1 — |(S''^'^)p, where S is the 
scattering matrix, S""^ = 6°^^ — iZ"^^, and Z'^^ is given 
by Eq. ([5]). Again, we take a = 1,2 to be the channels 
corresponding to the left and right scattering states, re- 
spectively. From Ref. [21| we know that in the random 
matrix framework and at the center of the energy spec- 
trum 



4k 



T1.2 



1,2 



(l + '«l,2)2' 



(27) 



where ki^2 is the effective coupling to channel 1 and 2, 
respectively. In the special case when the internal system 
is described by a GOE random matrix. 



Kl,2 



7r7i,2 

2A^i:> 



(28) 



where D is the mean level spacing at the center of the 
spectrum and N is the dimension of the internal system. 
As before, in the presence of asymmetry, 71 = 7 and 

72 = 7/9- 

For asymmetric coupling, maximum transmission is 
achieved when ti ~ T2 [18]. At this special point, the 
system with asymmetric coupling behaves as a system 
with symmetric coupling. This happens in the trivial 
case Ki = K2 (where the coupling is symmetric to be- 
gin with) or for ki = 1/k2- Since in the Id case we 
have ND « 27rr2 at the center of the energy band and 
for moderate disorder, we find from Eqs. (P7)) and (^5)) 
that the maximum transmission should be achieved when 
7/ri = 'i^/q, which is off by a factor of two from the value 
given by Eq. ((TU]). 

In order to understand the origin of this difference, we 
rederive below Eqs. ((27)) and (|28|) in a slightly different 
way, following the approach of Refs. [2l|. The scattering 
matrix averaged over the ensemble of random realizations 
is given by 



(S) 



1 - i{K) 



1 



(29) 



to leading order in 1/iV [H, [2l|. Here S"^ and K"'' are 
matrices in the channel space, with the explicit expres- 
sion for K'^'' given below in Eq. ^U^. Eq. (gH]) is vahd, 
under the assumption that the internal system can be de- 
scribed by random matrix theory, when the elements of 
iiT-matrices in the numerator and the denominator of Eq. 
(P5)l are effectively uncorrelated since their correlations 
lead only to corrections no larger than ~ 1/iV. In the 
Id Anderson model, the above assumption breaks down 
both in the limit of very weak disorder, where the internal 
system approaches integrability, and also for extremely 
strong disorder, where we enter the localized regime. 
The X-matrix in channel space can be written as 



K^^iE) 



Ij: 



E — En 



(30) 



where real energies En are the eigenvalues of the closed 
system. In our case of two channels, B^^ — ^J^{l\ipn) 
and Bn — y/j2{N\ipn) are the transition amplitudes for 
the eigenstate ]-0„) to the left and right leads, respec- 
tively. The sum in Eq. (pO)) contains the eigenvalues of 
the resolvent 1/{E— H) of the closed system. For energy 
E inside the spectrum of H we should understand it as 
a limiting value, E ^^ E + ii). Using the identity 



1 



E - En 



iO 



= P.v. 



1 



E — En 



i7rS{E~En), (31) 



and replacing the summation with the integral, we can 
compute the average .K-matrix. The principal value part 
is a smooth function of energy that vanishes in the middle 
of the spectrum; as a result, in this vicinity 



{K'') 



\Bi 



__^\^)p{E^O)^-tKu 



(32) 



where p{E = 0) = 1/D is the density of states at the 
center of the spectrum. Thus we have 



ttIB 



«;i,2 = 



B„=Ol 



2D 



(33) 



which we can use in Eq. (j27l) . In particular, if the eigen- 
state components are assumed to obey random matrix 
statistics, we have \B^ ^qP — 712/^, and we recover 
Eq. (Eg). 

We analyze the statistics of the Id Anderson model 
eigenstates in Fig. |8l where we plot the ensemble- 
averaged value of |(l|V'n)P (the probability overlap of 
the eigenstate IV'n) with a state localized at the left edge 
of the chain, that for a weakly open system would de- 
termine the width distribution for intrinsic states) as a 
function of En for different strengths of disorder. Evi- 
dently, the components of the eigenstates do not follow 
random matrix theory for moderate disorder. Near the 
center of the energy band we have (|(l|V'n)P) ~ 2/iV for 
W/n = 0.5 and W/n = 1. Then from Eq. dM]) we obtain 



Kl,2 



71.2 

217 



(34) 



and 7cr/^ ~ 2^^, in agreement with our findings in the 
previous sections. The values of 7/ri at which we have 
perfect tunneling probability, 7/17 = 2, where ti = 1, 
and 7/17 = 2(j, where T2 = 1, coincide precisely with the 
values of 7 at which the two super-radiant transitions 
occur, see the discussion in Sec. lIVI The fact that perfect 
tunneling probability r is reached at the super-radiant 
transition has been pointed out in Refs. [21J . 
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FIG. 8: (Color online) Overlap probability of an eigenstate at 
energy En with the edge 1 1) of the chain in the Id Anderson 
model for A'^ = 100. The result is averaged over 10* random 
realizations and plotted versus energy for different disorder 
values W (in units f2 = 1). The theoretical value for a GOE 
random matrix is indicated by a dashed horizontal line. 

From Fig. [5] we see that for large values of W we have 
(|(l|'0n)P) ~ 1/-^: but in this regime the eigenstates are 



strongly localized and we no longer expect Eq. (|29|) to be 
valid. The dip in the value of (|(I|?A„)p) at the center of 
the energy band shown in Fig. [8] is consistent with the 
analysis of Ref. [17], where it was pointed out that the 
localization length is shorter at the center of the energy 
band than for surrounding energies. 
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FIG. 9: (Color online) Left and right tunneling probabilities 
ri,2 are plotted vs 7/f2 for three different disorder strengths 
W, in units SI = 1, for the system size A'^ = 100 and asym- 
metry parameter q = 4. The numerical results (symbols) are 
obtained by averaging over 10'^ disorder realizations and over 
an energy interval —0.1 < E < 0.1 around the center of the 
spectrum. These are compared with the theoretical results 
(solid curves) obtained from Eqs. (|34p and H27|) . The critical 
value of 7/12, Eq. ([lOp . is indicated by a vertical solid line 
in each panel. In the upper and middle panels, the values at 
which perfect transmission is reached, 7/12 = 2 and j/il = 2g, 
are shown as vertical dot-dashed lines. 

In Fig.|9]we compare our numerical results for the tun- 
neling probability r with Eqs. ((77)) and (|34l) . showing 
reasonable agreement for moderate disorder W/il = 0.5, 
1. For W/H. = 10, clear deviations from the analytical 
expressions are visible, due to the fact that the assump- 
tions of random matrix theory break down at very strong 
disorder. 

These results show that a blind application of random 
matrix results would lead to incorrect conclusions for the 
Id Anderson model. Our empirical expression for the 
tunneling probability works well for moderate disorder. 
In this regime, we regain the critical value of the cou- 
pling strength (jlOp for which maximum transmission is 
achieved. 



VI. MULTI-DIMENSIONAL CASE 

In this Section, we discuss the higher-dimensional 
cases. Only selected results will be shown, sufficient 
to demonstrate that the general behavior of the maxi- 
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mum transmission found in Id systems can be extended 
to higher dimensions. 

The open model in dimensions greater than one con- 
sists of an array of sites with associated energy levels. 
Neighboring sites are coupled to each other by the tun- 
neling coupling J7, as in the Id case. In higher dimen- 
sions, we have many ways to couple the system to exter- 
nal leads. In the 2d case, where the number of sites in a 
rectangular array is N = AI x L, we couple each of the 
M sites on the left side to a separate left lead and each 
of M sites on the right side to a separate right lead, see 
Fig. [TUl The coupling amplitude for each of the M left 
leads is 71 = 7, while the coupling amplitude to each of 
the M right leads is 72 = 7/g. Similarly, in the 3d case, 
with N = M X M xL sites, each site on an M x M face is 
coupled to a lead. In this geometry, each lead represents 
a channel. Such an open model can describe a variety of 
physical systems, such as an array of quantum dots, or a 
particle trapped in a lattice potential. 
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FIG. 10: Two-dimensional open model used in this paper, 
with N = M X L sites coupled to M incoming left channels 
and M outgoing right channels. Specifically, in this Figure 
we show M X L = 5 X 5 sites coupled to 2M = 10 channels. 
At each site there is a bound state with energy Eo coupled to 
its nearest-neighbor sites through a tunneling amplitude fl. 

The diagonal part of the effective Hamiltonian for this 
system can be written as T-Lu = Eq + SEq for sites i that 
are not coupled to leads, and Hu = Eq + SEq + (5i,2 ~ 
(i/2)7i^2 for sites coupled to the left or right leads, re- 
spectively. As in the Id case, we set the center Eq of 
the energy band to zero; SEq is a random variable uni- 
formly distributed in [-~W/2,+W/2], and PF is a disor- 
der parameter. As before, ^1,2 and 71.2 represent the 
energy shift and decay probability (inverse lifetime), re- 
spectively, induced by the coupling to the left and right 
leads. In the following we again neglect the energy shift, 
so that Si^2 = 0. Finally, for the off-diagonal matrix ele- 
ments, i j^ j, we have Hij — Hji = VL \i the sites i and j 
are neighboring and T-Lij — otherwise. 



Depending on the degree of disorder, different trans- 
port regimes are possible in the multidimensional case. 
The ballistic regime is defined by the condition L < I, 
where L is the system length and I is the mean free 
path. The diffusive regime is determined by the condi- 
tion I < L < Lioc, where as before Lioc is the localization 
length. Finally, for Lioc < L, we are in the localized 
regime. The mean free path and the localization length 
both depend strongly on the energy interval under con- 
sideration and on the disorder strength, see the discussion 
in Refs. :22, 23]. 

In our multidimensional arrangement, a double super- 
radiant transition again occurs as in Id, which we do not 
discuss in detail here. In this Section we will focus on 
the dependence of the maximum conductance on the cou- 
pling strength to the leads. The dimensionless conduc- 
tance G, which is proportional to the total transmission, 
can be computed using the Landauer formula [I^, [23] , 
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2J\/ 
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(35) 



a=l b=M+l 



Here Z°-^{E) is the transmission amplitude between chan- 
nels a and b, see Eq. ([5]), that can be computed from the 
effective Hamiltonian. 

We will not consider here the case of a small pertur- 
bation to an integrable Hamiltonian, which displays very 
system-specific behavior, and start by analyzing the case 
of moderate disorder. In Fig. [TTJ we illustrate the be- 
havior of the average conductance in the diffusive regime 
for a quasi-ld lattice of M x L = 10 x 100 sites with 
W/Q = ^3/4, upper panel, and for a square 2d lat- 
tice of 30 X 30 sites with VF/fi — 2, lower panel. In 
both cases, the maximum conductance is obtained at the 
critical value of the coupling given by Eq. (flUl) . In the 
upper panel, we compare our numerical results with an- 
alytical expressions obtained in the context of random 
matrix theory 18*. The numerically computed left and 
right transmission coefficients as functions of ^/fl were 
used to evaluate the analytical expression for the aver- 
age conductance given in Ref. jlSl ]. In the upper panel 
of Fig. Ill[ the numerical results (symbols) are in good 
agreement with the analytical results (solid curves) for 
a quasi-ld system in the diffusive regime. To the best 
of our knowledge, there are no known analytical results 
for the 2d case in the diffusive regime with asymmetric 
coupling (lower panel). 

At the critical value of j/fl given by Ec^. (|10p . each 
conductance curve for the asymmetric case intersects the 
conductance curve for the symmetric case {q = 1). This 
is similar to what we have observed in the Id case (see 
Fig. [3]), and is also consistent with the results obtained 
in the context of random matrix theory. As discussed 
in Sec. |Vl the maximum conductance is achieved when 
Ti — T2. In this case a system with asymmetric coupling 
behaves as a system with symmetric coupling. 

In the regime of very strong disorder, the critical value 
at which we have maximum conductance becomes depen- 
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again given by Eq. ((TO)l . A detailed comparison between 
random matrix results and the Anderson model will be 
presented elsewhere. 
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FIG. 11: (Color online) The average dimensionless conduc- 
tance as a function of the coupling strength 7/r2. In the upper 
panel we consider the quasi- Id case of 10 x 100 sites in the dif- 
fusive regime, with disorder W/Q, = -^3/4, so that L/l ^ 3. 
The analytical results of random matrix theory [13] (solid 
lines) are compared with our numerical results (symbols). In 
each case, an overall multiplicative factor in the analytical 
expression J_0| has been obtained by fitting to the numeri- 
cal data. In the lower panel, we consider the 2d case with 
30 X 30 sites and W/Q = 2, so that L/l « 4. Here symbols 
refer to our numerical results and the dashed lines simply 
connect the symbols (random matrix results are not available 
for asymmetric 2d systems in the diffusive regime). Vertical 
dashed lines indicate the critical values 7/f2 = 2y^, where 
the maximum conductance is expected for each asymmetry 
value q. All numerical results are obtained by averaging over 
200 realizations of disorder and over 100 different energies in 
the interval —0.1 < E < 0.1 around the center of the energy 
band. 



dent on the mean level spacing, D. Similarly to the Id 
case we have 



7cr oc ND . 



(36) 



Since in this regime D sa W/N, we find that the critical 
coupling for the maximum conductance is proportional 
to the disorder strength W, see Fig.[T2] 

For strong disorder, we enter the localized trans- 
port regime, where the transmission is log-normally dis- 
tributed [l^, so that, similarly to Fig. |4l it is more con- 
venient to consider the average of In G rather than the 
average of G. In Fig. [12] we plot the average logarithm 
of the conductance as a function of 7/^2 in 2d, upper 
panel, and in 3d, lower panel. As indicated by the ar- 
rows, the critical value of 7/ri at which the conductance 
is maximized is in both cases proportional to the mean 
level spacing, which in turn is proportional (for strong 
disorder) to the disorder strength. For moderate disor- 
der in either 2d or 3d, see e.g. W — 2 in the upper 
panel and M^ = 5 in the lower panel, the critical value 
of 7/il for which we have the maximum conductance is 




FIG. 12: (Color online) Average logarithm of the conduc- 
tance as a function of the coupling strength y/Q for the mul- 
tidimensional case, with SI = 1, asymmetry parameter g = 4, 
and different disorder strengths W . For the sake of clarity, 
we plot (InG) — (InG)o so that all data start from the same 
point, where (InG)o is defined as in Fig. [l] In the upper 
panel, we show the 2d case with 20 x 20 sites, while in the 
lower panel we show the 3d case, with 8x8x8 sites. The 
dashed vertical line shows the critical value of the coupling, 
Eq. (|10|l . at which we have maximum conductance for mod- 
erate disorder. The arrows indicate "f/^l = 0.76ND/Q for the 
2d case, upper panel, and 0.56ND/fl for the 3d case, lower 
panel. The numerical factors 0.76 and 0.56 have been found 
from fitting. All data have been obtained by averaging the 
conductance over 10* ensemble realizations at the center of 
the energy band, E = 0. 



VII. CONCLUSION 

We have analyzed coherent quantum transport 
through a finite sequence of potential barriers with asym- 
metric coupling to the external environment. As the 
coupling to the environment is varied, transport prop- 
erties are greatly affected. In a previous paper j4i], the 
super-radiant transition that occurs in this paradigmatic 
model, at a critical value of the coupling, was studied for 
the case of symmetric coupling to environment. Here, 
with the aid of the effective non-Hermitian Hamiltonian, 
we show that for asymmetric coupling a double super- 
radiant transition occurs, as compared with a single tran- 
sition in the symmetric case. 

The super-radiant transitions have important conse- 
quences for the observable resonance structure. In par- 
ticular, the number of resonances decreases by one after 
each transition. Focusing on the behavior of transmission 
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near the energy band center, see Fig. [71 we demonstrate 
that a sharp change in the structure of the resonance oc- 
curs in correspondence with the two super-radiant tran- 
sitions. This change can be characterized by the behav- 
ior of the transmission maxima and minima, which we 
describe analyticahy as functions of the couphng to the 
environment and of the asymmetry of this couphng. As 
far as we know, these features of the structure of reso- 
nances as a function of the couphng strength to the leads 
have not previously been reported in the literature. 



Appendix A: Analytical Derivation of the 
Transmission in the Ordered Case 



To derive closed analytical expressions for the sums 
(13), it is convenient to double them, making the phase 
Lpn = WK / [N + 1) run over the whole circle. Then 

1 1 2Ar-|-l .2/ ^ 



Maximum transmission through the system is reached 
at a coupling ^/Vt — 2y/q, where q is the asymmetry pa- 
rameter. This coupling is equal to the geometric mean 
of the coupling strengths associated with the two tran- 
sitions. We show that this result does not follow from 
random matrix theory, as usually assumed in statistical 
theories of quantum transport. Moreover, for very strong 
disorder we show that the coupling at which transmission 
is maximized is proportional to the disorder strength. We 
also find that the latter results remain valid in higher 
dimensions. Specifically, we analyze the average con- 
ductance as a function of the degree of opening, of the 
asymmetry, and of the strength of disorder in the multi- 
dimensional cases: quasi-ld, 2d, and 2>d. In the quasi-ld 
case, we compare our results with analytical expressions 
obtained in the context of random matrix theory. We 
demonstrate the validity of our results in 2d and 3d, in 
both diffusive and localized transport regimes, where, to 
the best of our knowledge, no analytical results as a func- 
tion of the asymmetry of the opening are available in the 
literature. 

The results presented here are based on an approach 
originally formulated in the framework of nuclear reac- 
tion theory. Now we understand that they reflect gen- 
eral properties of quantum signal transmission. There- 
fore they might be of relevance for numerous applica- 
tions. The sequence of potential barriers is a paradig- 
matic model for coherent quantum transport. A better 
understanding of this transport regime is essential for 
the development of information technology using differ- 
ent nanoscale systems with complex geometry, including 
quantum dots, photonic crystals, and molecular wires. 
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where e — E/2Q.. Now the sum over n spans all the roots 
Zn = exp(i(/?„) of the equation 2;2Ar+2 _ i — q, so that 

*^ ^ 4f7 2iV + 2 ^/ ^ z„(z2_2ez„ + l) 

n— 

(A2) 
Splitting the summand into simple fractions and using 
the identities 
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we can write down the original sums in closed form. 
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(A7) 
where z± — e ± \J e^ — 1 are the roots of the quadratic 
equation z^ — 2ez -1-1 = 0, and z+z_ = 1. In the 
asymptotics of large N and energy outside the Bloch 
band, |e| > 1, we have either \z^\ ^ 1 ^ \z^\ or 
|z:^| ^ 1 :^ |z:^|, for e > 1 or e < — 1 respectively, while 
inside the band we have z± = exp(±i/3). This leads to 
Eqs. (HZl) - dH]) in the main text. 

To calculate the transmission amplitude (J12p or trans- 
mission probability (|15|) . we also need 



Pl{e)-Pl{,) 



1 sin[(Ar - l)/3] 
4fi2 sin[(Ar + l)/3] ' 



(A8) 



valid for all energies inside the band, |e| < 1. For A^ ^ 1, 
we can consider averaging over a small energy interval, 
similarly to Eq. (P5)) . where e = cos/3 is constant while 
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the trigonometric functions depending on Nfi change 
from —1 to +1. This can be done formahy with the 
substitution /3 — ^ /3 + i/N , where ^ is of order one. Then 
the average transmission T = T^"^ = jZ^^p is given by 
the integral 



yi 



4?7iry2 



d{NP) 



(A9) 



TV $(iV/3) ' 
where 771^2 = 7i. 2/2^2 and 

$ (y) = A^ sin^ y + B^ cos^ y + 2C sin y cos y , ( AlO) 



A^ = (r/i + r]2y + (1 - mm) cos" /3 



B^ = (l+?7i7?2)"sin"/3 



C = (1 — 7/^772 ) cos /3 sin /3 . 



(All) 



(A12) 



(A13) 



After changing the integration variable to a; = tan(A^/3), 
the integral becomes 



Ii(3) = - 

TT 



dx 



1 



A^x^ + 2Cx + 52 V^2^2 _ (^2 ■ 

(A14) 



This leads to the final result for the average transmission 
at any energy inside the band, 



7fi2_ 4?7i772sin^ 



2^12 _ 



im +'72)(i + riim) 



(A15) 



where sin/3 = \/T~~?, to be compared to Eq. ((23|) . By 
means of slightly more complicated integrals, we derive 
[here cos(2;3) = 2e^ - 1] 

yTT^4 2 ??i(l + '7i'72 + ??|) + m cos(2/3) + (1 + ?7|) sin/3 
''^ (?7i+ry2)(l + r7i»72)(l + '7? + 277isin/3) 

_ (A16) 

and similarly for T^^. This, along with Eq. (fT2|) . allows 
one to check the unitarity condition. 



rj.11 ^rj.12 ^ _2ImZ" 



(A17) 



which remains valid after averaging, where both sides are 
equal to 



4?7i (?7i + sin /3) 
1 + r]f + 2rji sin ^ 



(A18) 
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